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Abstract 

Little is known about the microevolutionary processes shaping within river popula- 
tion genetic structure of aquatic organisms characterized by high levels of homing 
and spawning site fidelity. Using a microsatellite panel, we observed complex and 
highly significant levels of intrariver population genetic substructure and Isolation- 
by-Distance, in the Atlantic salmon stock of a large river system. Two evolutionary 
models have been considered explaining mechanisms promoting genetic substruc- 
turing in Atlantic salmon, the member-vagrant and metapopulation models. We 
show that both models can be simultaneously used to explain patterns and levels of 
population structuring within the Foyle system. We show that anthropogenic factors 
have had a large influence on contemporary population structure observed. In an 
analytical development, we found that the frequently used estimator of genetic dif- 
ferentiation, f sTi routinely underestimated genetic differentiation by a factor three 
to four compared to the equivalent statistic Jost's D^st (lost 2008). These statistics 
also showed a near-perfect correlation. Despite ongoing discussions regarding the 
usefulness of "adjusted" f st statistics, we argue that these could be useful to iden- 
tify and quantify qualitative differences between populations, which are important 
from management and conservation perspectives as an indicator of existence of bi- 
ologically significant variation among tributary populations or a warning of critical 
environmental damage. 



doi: 10.1002/ece3.32 



Introduction 

Over the past decade, microsatellites have been used widely 
as the molecular marker of choice for studies of popula- 
tion structure in Atlantic salmon {Salmo salar). The spatial 
scales of these studies have varied from intercontinental (King 



et al. 2001), intracontinental (Spindle et al. 2003), to regional 
(e.g., Fontaine et al. 1997; DiUane et al. 2007). Given their 
greater discriminatory power in comparison to other molec- 
ular markers (e.g., protein electrophoresis, Heggberget et al. 
1986), microsatellites are particularly effective for studies ex- 
amining population genetic structure on small spatial scales 
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Table 1. Comparison of sampling design and geography of studied area in recent studies examining evolutionary models explaining within-river 
population structure in Atlantic salmon {Salmo salar). 





Garant et al. 


Primmer et al. 


Vaha et al. 


Dillane et al. 






(2000) 


(2006) 


(2007) 


(2008) 


This study 


System studied 


Sainte-Marguerite 


Varzuga 


Teno/Tana 


Moy 


Foyle 


Number of loci 


5 


17 


32 


12 


7 


Mm m hoc nt <; a m n 1 o c 


683 


392 


792 


1 606 


1 086 


Geomorpliological feature 


Linear system 


Linear system 


Linear system 


Lacustrine system 


Dendritic system 


Catchment size 


llUkm^ 


9510 km^ 


16,386 km^ 


2000 km^ 


4450 km^ 


Life stage sampled 


Fry 


Parr 


Adult 


Fry/Parr 


Fry/Parr 


Temporal genetic stability test 


Yes 


No 


No 


Yes 


Yes 


Evolutionary model 


Meta population/ 


Member Vagrant 


Metapopulation 


Metapopulation 


Metapopulation/ 




Member Vagrant 








Member Vagrant 


Climatic zone 


Subarctic 


Arctic/Subarctic 


Arctic/Subarctic 


Cool temperate 


Cool temperate 


Anthropogenic disturbance 


Moderate (hydro power) 


Limited 


Limited 


Substantial (drainage) 


Substantial (drainage) 



(e.g., Wright and Bentzen 1995; O'Connell and Wright 1997). 
A small number of studies (Garant et al. 2000; Primmer et al. 
2006; Vaha et al. 2007; Dillane et al. 2008) have focused 
on within-river population structure in Atlantic saknon, and 
these differed considerably, for instance, in the number of mi- 
crosatellite markers used, total number of samples, temporal 
replicate sampling, river catchment size, geomorphological 
structure, and life-history stage targeted (Table 1). 

A common finding of these studies, to a greater or lesser 
degree, was the existence of within river genetic population 
structure. For instance, in the west of Ireland, within-river 
genetic variation has been reported in the River Moy. Dillane 
et al. (2008) observed a considerable degree of genetic vari- 
ation vnthin and between tributaries of this comparatively 
large and complex river system, with the presence of two large 
lakes as a significant landscape feature, explaining the popula- 
tion structure observed. Over the past few years, a number of 
similar studies have attempted to identify the best underlying 
evolutionary models responsible for observed patterns of ge- 
netic structuring of wild Atlantic salmon populations (Garant 
et al. 2000; Primmer et al. 2006). Two evolutionary models 
have often been proposed to describe the underlying cause 
of such structuring; the "metapopulation" and the "mem- 
ber/vagrant evolutionary models." Under the metapopula- 
tion model, local subpopulations are linked by gene flow and 
experience recurring extinction-recolonization events due 
to environmental instability (Hanski 1998). Consequently, 
strong genetic structuring and locally adapted gene pools are 
prevented from being formed. A temporally stable genetic 
population structure and a positive correlation between ge- 
netic differentiation and geographical distance (Isolation-by- 
Distance, IBD) are unlikely to occur in a true metapopulation. 

The member-vagrant evolutionary model predicts that 
precise homing restricts gene flow, thus favoring the forma- 
tion of locally adapted gene pools and subsequent strong and 
temporally stable genetic structuring among populations. 



with a clear IBD signal (e.g., Sinclair 1988; Garant et al. 2000). 
Primmer et al. (2006) reported on strong, weU-supported 
population structuring in combination with clear IBD in the 
Atlantic salmon population of the Varzuga River (Russia) and 
concluded that this was in accordance with the predictions of 
the member-vagrant evolutionary model though, one of the 
specific predictions of the member-vagrant model, namely 
temporal stability, was not tested by Primmer et al. (2006). 
Garant et al. (2000) reported weak genetic structuring and 
lack of IBD in the Canadian Sainte-Marguerite River, and 
found evidence for temporal stable genetic substructuring 
at four of seven sampling locations, concluding that both 
the member-vagrant and the metapopulation models can be 
used to describe the genetic structure observed. 

Knowledge of the evolutionary model most suitable to 
explain observed genetic differentiation in salmonid popula- 
tions can be of great benefit in designing and implementing 
effective management and conservation strategies to ensure 
long-term sustainability of fish stocks (e.g., to identify and 
assess the biological and genetic health status of individual 
population units). Given the lack of consensus as to which 
evolutionary model best explains patterns of genetic differen- 
tiation in wild Atlantic salmon populations, and that studies 
focusing on within-river genetic structure in this species re- 
main relatively scarce, more studies are warranted. 

The River Foyle system, in the northwest of Ireland, was, 
until recently, one of the most productive wild Atlantic 
salmon fisheries in Western Europe. Annual commercial 
catches of salmon in the estuary and the adjacent sea area have 
ranged between 12,176 and 44,000 salmon between 1997 and 
2006 (Loughs Agency Annual Report and Accounts 2007), 
while annual angling catches for the River Foyle were esti- 
mated to be in the region of 10,000 fish. There has been a ces- 
sation of commercial salmon fishing, seaward of Lough Foyle 
since 2007 in order to reduce the potential for the exploitation 
of mixed stocks outside of the Foyle, thus changing the focus 
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of exploitation toward a single river fishery with an identifi- 
able surplus. As a result, commercial catch figures dropped 
to 12,176 fish in 2006 and 5372 in 2007 (Loughs Agency 
Annual Report and Accomits 2007). From 2010, commercial 
fishing on the River Foyle has been completely banned as a 
conservation measure. 

The Foyle system, in contrast to the Moy, Varzuga, Teno, 
and Sainte-Marguerite (Table 1), is a temperate (geograph- 
ically), medium-sized, geomorphologically complex system, 
comprising an extensive dendritic network of major and 
minor tributaries rivers with only one significant lake (Lough 
Derg), which is a headwater lake upstream of available salmon 
spawning habitat and thus unlikely to contribute to popula- 
tion structure, representing a variety of habitat types and 
flow regimes. Returning adult spawners in the River Foyle 
system consist mainly of grilse (one-sea winter; ISW) with 
one of the larger tributaries (River Finn) having a major com- 
ponent of spring-running multi-sea-winter (MSW) salmon. 
Stocks in several of the tributaries display characteristically 
late running grilse components. Thus, this complex river sys- 
tem provides an excellent opportunity to examine further 
evidence for the influence of microevolutionary processes 
shaping contemporary population genetic structure, which 
both contrast and complement previous studies. 

The specific aims of this study are to: ( 1 ) determine spatial 
genetic structure in the Foyle system, (2) asses which evo- 
lutionary model (member-vagrant or metapopulation) best 
fits the observed results, and (3) compare and discuss results 
with other such studies from the recent literature. 



Materials and Methods 
study site 

The Foyle River system (55°00'N; 07°20'W) is located in the 

northwestern corner of the island of Ireland, forming part of 
the land border between Northern Ireland (United Kingdom) 
and the Republic of Ireland (Fig. 1). The Foyle River system, 
together with the Faughan River, falls within a single man- 
agement jurisdiction, the cross-border Loughs Agency. The 
river has a catchment area of 4450 km^. The main tributaries 
of the Foyle are the rivers Finn and Derg to the west, and the 
OwenkUlew River in the eastern part of the catchment. The 
Foyle main stem is known under different names throughout 
its course. Downstream of the confluence of the Drumragh 
River and the Camowen River, the main stem is called River 
Strule until it meets the River Derg downstream of which the 
main stem is known as River Mourne. After the River Finn 
joins the River Mourne, the main stem is called River Foyle 
until it drains into Lough Foyle. The Faughan River does not 
join the Foyle main stem, but directly enters Lough Foyle to 
the east of the main river. 

Sampling strategy 

Sampling was conducted throughout the Foyle system, in- 
cluding the adjacent River Faughan, with an average water- 
way distance between the sampling locations of 72 ± 28 (SD) 
km and a range of 17.8-132.4 km. Following consultation 
with local fishery managers, sample coverage included most 
known spawning and nursery areas in all major tributaries. 




Figure 1. Map of Ireland (inset) and the 
Foyle River system. Symbols indicate 
location of sampling sites. Abbreviations are 
as given in Table 2. 
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Table 2. Details of sampling sites in the Foyle River system. 



Site 


Abbreviation 


Sample year 


Coordinates 


Region 


N 


Faughan 


Fau 


2006 


54°57'N; 07M7'W 


Northeast 


96 


Faughan 


Fau_99 


1999 


54°57'N; 07°17'W 


Northeast 


45 


Deele 


Dee 


2007 


54°51'N; 07°41'W 


Northeast 


48 


Deele 


Dee_99 


1999 


54°51'N; 07°41'W 


Northeast 


33 


Burn Dennet 


Den 


2006 


54°B3'N; 07°23'W 


Northeast 


96 


Struie (Cappagh Burn) 


Str 


2007 


54°40'N; 07°18'W 


Northeast 


48 


Owenreagh East 


Owe 


2006 


54°41'N; 07"05'W 


Northeast 


48 






2007 




IMUl LI icci^L 


48 


Finn 


Fin 


2006 


54°47'N; 07°45'W 


West 


96 


Reelan 


Ree 


2006 


54°49'N; 08°03'W 


West 


96 


Fairy Water 


Fai 


2007 


54'=38'N; 07°23'W 


West 


48 


Derg 


Der 


2006 


54°39'N; 07°44'W 


West 


96 


Owenreagln South 


Ows 


2006 


54°32'N; 07°29'W 


Southeast 


96 


Glenelly 


Gle 


2006 


54°46'N; 07°13'W 


Southeast 


96 


Drumragh 


Dru 


2007 


54°28'N; 07°17'W 


Southeast 


96 



Sampling by electrofishing conducted during 2006 and 2007 
captured juvenile Atlantic salmon from multiple year classes; 
48 or 96 fish per sampling location with ideally a 50/50 ratio 
of 0+ (fi-y) and 1+ (parr) fish. A total of 13 sampling sites 
were visited. Among these, two temporal replicate samples 
originally taken at two sampling sites (Deele and Faughan 
Rivers) in 1999 were included to test for temporal variation. 
In total, 1086 individual tissue samples were collected for this 
study (Fig. 1; Table 2). The idealized sampling ratio was not 
achieved at every sampling location, indeed at some sam- 
pling sites visited during 2007 no 0+ sample was taken due 
to the fact that at the time of year of the visit, the O""" Atlantic 
salmon could not yet be distinguished from 0"*" brown trout 
(S. trutta). The aim was to collect 96 samples per sampling 
location. At some of the (smaller) tributaries, this number 
was not reached due to smaller numbers of juvenile salmon 
present. In these cases, the number of samples taken was 
reduced to 48. 

The sampling of multiple year-classes combined with tak- 
ing the samples from river stretches of several hundred me- 
ters minimized the likelihood for the "AUendorf- Phelps ef- 
fect" (Allendorf and Phelps 1981), whereby false statistically 
significant genetic differentiation among sampling sites can 
result from the sampHng of the offspring of a very limited 
number of individuals that do not represent the full extent of 
genetic variation present in the population. 

Tissue samples from the juvenile Atlantic salmon were 
taken by removing the adipose fin before releasing the fish 
in the case of the parr samples, or by sacrificing the whole 
fish for the fry samples. In both cases, the tissue sam- 
ples were preserved in molecular grade 99% ethanol until 
analysis. 



Genetic analysis 

DNA was extracted by using a ChelexlOO (Bio-Rad, Hemel 
Hempstead, United Kingdom) protocol (Estoup et al. 1996) 
on small amounts ( ~ 3 mm^ ) of fin tissue. Extracted DNA was 
kept frozen at -20° C until amplification using a polymerase 
chain reaction (PGR). Samples were screened for a suite 
of seven highly polymorphic microsateUite loci; Ssspl605, 
Sssp2201, Sssp2210, Sssp2216, SsspG7 (Paterson et al. 2004), 
Ssal97, andSsa202 (O'Reilly etal. 1996). The forward primer 
of each primer pair was labeled at the 5' end with a fluores- 
cent dye: 6FAM (blue) for Ssal97, Ssspl605, Sssp2201, and 
Sssp2210 loci, VIC (green) for loci Sssp2216 and SsspG7, 
NED (yellow) for locus Ssa202. 

The individual PGR amplification reactions for loci 
SsspG7, Ssa202, and Ssa 197 consisted of 1 /il of extracted tem- 
plate DNA, 0.2 mM dNTP's (Invitrogen, Paisley, Scotland, 
United Kingdom), 0.7 x reaction buffer (13 mM Tris-HGl 
[pH 8.4], 33 mM KGl; Invitrogen), 1.5 mM MgGlj, 0.3 /^M 
forward and reverse primer (0.15 jxM in the case of Ssal97), 
0.13 U Taq DNA polymerase (Invitrogen), and ultraPURE 
H2O (Gibco, Paisley, Scotland, United Kingdom) to a total 
reaction volume of 1 5 /xl. For the remaining four loci, the PGR 
reaction consisted of 1 /xl of extracted template DNA, 0.2 mM 
dNTP's (Invitrogen), Ix reaction buffer (20 mM Tris-HGl 
[pH 8.4], 50 mM KGl; Invitrogen), 5 mM MgGl2, 0.38 /xM 
of the forward and reverse primer for Ssspl605, 0.95 jOtM for 
Sssp2201, 0.17 /xM for Sssp2210, 0.88 /xM for Sssp2216, 0.1 U 
of Taq DNA polymerase, and ultraPURE H2O to a total reac- 
tion volume of 10 /xl. The PGR, in which Sssp2210/Sssp2201 
and SsspG7/Ssal97 were multiplexed, was performed using 
the following amplification cycles: denaturing at 95°C/2 min 
followed by 33 cycles of 95°G/45 sec, 58°G/45 sec, 72°C/45 sec. 
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which was foUowed by a final cycle of 95°C/45 sec, 58°C/45 
sec, 72°C/5 min on a Thermo Hybaid PGR Express thermo- 
cycler (Thermo Hybaid, Franklin, MA, United States). 

The PGR products were prepared for electrophoresis on 
an ABI 3130 Automatic DNA Analyzer (Applied Biosystems, 
Paisley, Scotland, United Kingdom) according to the man- 
ufacturer's instructions using ABI GS-500 LIZ internal size 
standard. Alleles were scored and sized using the automated 
bin function of the ABI GeneMapper v4.0 software, but were 
manually checked and verified. 

Quality control 

Microchecker v2.2.3 (van Oosterhout et al. 2004) was used to 
test for the possible presence of nuU-alleles and allele scoring 
errors. To further check for geno typing errors, 118 individuals 
( 1 1% of the total data) from two randomly chosen sampling 
locations were genotyped again following the same protocols 
as in the original analysis. The genotypes were compared to 
the previous genotypes and allelic mismatches were counted 
to establish an error rate for the genotyping in this study. 

Statistical analysis 

Genetic diversity in all samples was calculated; num- 
ber of alleles (ATa), allelic frequencies, observed (Hq), 
and expected (H^) heterozygosity, using the program 
GENEPOP v3.4 (Raymond and Rousset 1995). Deviations 
from Hardy-Weinberg equilibrium (HWE) for each sam- 
ple (Guo and Thompson 1992) were calculated using the 
Genepop 3.4 software (Raymond and Rousset 1995). 

Temporal genetic stability, between cohorts (i.e., short- 
term temporal stability) and between samples Deele and 
Faughan and their 1999 equivalents, was assessed by testing 
for pairwise heterogeneity in allelic frequency distributions 
between temporal samples with GENEPOP v3.4. The same 
program was also used to test for genotypic linkage disequi- 
librium across putative populations, to confirm independent 
allelic segregation among different marker loci. 

Genetic differentiation between samples was quantified 
by calculating pairwise f st values, or variance in allele fre- 
quencies (Weir and Cockerham 1984), using FSTAT v2.9.3. 
(Goudet 2001). The statistical significance of estimated val- 
ues was tested by permutation. The same programme, FSTAT 
v2.9.3., was also used to calculate allelic richness (Ar), which 
is the allele number per population corrected for differences 
in sample size. Corrections to account for multiple tests, se- 
quential Bonferroni correction, (Rice 1989) were carried out 
for all tests mentioned above. 

Genetic differentiation between samples was also calcu- 
lated using an unbiased estimator of Jost's D (Jost 2008) as 
implemented in the programme SMOGD (Crawford 2010); 
Dest. Recently, the validity of using Pst and related measures 
of genetic differentiation in population genetic studies based 



on microsateUite data has been questioned (e.g., Hedrick 
2005; Jost 2008; Whitlock 2011). The rationale is that Fst 
tends to be biased toward lower values, particularly in sit- 
uations where within-locus genetic diversity is high, as is 
often the case for highly polymorphic microsateUite data (see 
Hedrick 2005 and Jost 2008). Here, we compare Dest direcdy 
to Fsv to assess whether "traditional" measures of population 
differentiation underestimate differentiation in the current 
dataset, to what degree, and what possible consequences this 
has for further analysis and conclusions. For this purpose, 
possible correlations between f st and D^st were compared by 
means of a Mantel test (Mantel 1967). 

Genetic distances between samples were estimated us- 
ing the Cavalli-Sforza and Edwards (1997) Chord distance 
measure (Dc) as implemented in the programme POPULA- 
TIONS vl.2.30beta (LangeUa 1999). This measure of genetic 
distance was selected in favor of others as Dq is drift based, 
relatively independent from a specific mode of mutation of al- 
leles and fluctuations in effective population size, and is most 
likely to provide the correct phylogeny in closely related pop- 
ulations (e.g., Takezaki and Nei 1996; Goldstein and Pollock 
1997). The distance matrix was used to construct a neighbor- 
joining (NJ) phylogram with 10,000 bootstrap replicates over 
loci to obtain confidence estimates of the branching pattern, 
drawn with the programme TREEVIEW (Page 1996). To test 
for genetic IBD, pairwise linearized Fst and Dest estimates, 
calculated between sample sites, were plotted against geo- 
graphical distance using a Mantel test (Mantel 1967) imple- 
mented in IBDWS V3.15 (Jensen et al. 2005). The signifi- 
cance of observed associations were evaluated using 10,000 
permutations. 

To further investigate the level of population genetic 
structuring within the River Foyle system, the programme 
STRUCTURE v2.3. 1 (Pritchard et al. 2000) was also used. The 
algorithm buUt within this programme allows for the deter- 
mination of the most likely number of populations (clusters) 
explaining the data without requiring "a priori" definition 
of populations. The programme was set to run using the ad- 
mixture model with correlated gene frequencies. As STRUC- 
TURE uncovers only the uppermost level of structure in a 
dataset, Pritchard and Wen (2004) suggest running STRUC- 
TURE subsequently on each identified cluster separately to 
reveal any underlying structure; in effect a "hierarchical struc- 
ture analysis." This approach was followed in this study. 

To infer the true number of clusters (K), we followed 
Evanno et al's. (2005) "ad-hoc" statistic (AK) based on the 
rate of change in the log probability between consecutive K 
values, which ranged from K=ltoK = 20 with a burn-in and 
Markov chain Monte Carlo length of 10,000 each (100,000 
burn-in and 200,000 Markov chain Monte Carlo repUcates 
when analyzing for hierarchical within-cluster structure) for 
20 runs per ^T-value. The modal value of the AK distribu- 
tion is assumed to represent the true number of clusters or 
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populations in the data. The reliability of the replicate anal- 
ysis mentioned above was enhanced by use of the CLUMPP 
1 . 1 .2 ( Jakobsson and Rosenberg 2007) software package. This 
computer program's "Greedy" algorithm was used to find the 
optimal alignment of the results from the 20 replicate cluster 
analyses. A second method used to find the most likely num- 
ber of clusters was taking the maximum value of the mean 
value of the log likelihood (L[i<r)] of the data as the most 
likely value of K (e.g., Zeisset and Beebee 2001; DiUane et al. 
2008). 

The results from the STRUCTURE analysis were summa- 
rized using the programme Distruct 1.1 (Rosenberg 2004), 
which generates a bar plot of the STRUCTURE output data 
where each individual is represented by vertical colored line. 
Different colors code for each different K representing the 
individual's membership coefficients to each K. A hierar- 
chical analysis of spatial genetic diversity was performed 
using the method of molecular variance (AMOVA) imple- 
mented with the software Arlequin 3.1.1 (Excoffier et al. 
2005). 

Barriers to gene flow that shape patterns of genetic dif- 
ferentiation between populations were identified using the 
software package Barrier 2.2 (Manni et al. 2004), which uti- 
lizes Monmonier's maximum difference algorithm to identify 
such barriers. The robustness of the computed barriers was 
assessed by analyzing 1000 resampled bootstrapped Fsi ma- 
trices with the software package. 

Results 
Quality control 

No evidence for either large allele dropout or scoring er- 
rors due to stuttering was detected by MicroChecker. In the 
regenotyping exercise (involving 118 samples for 7 loci), only 
five genotyping errors were found; a genotyping error rate of 
0.3%. Three of these errors were allelic dropouts, the other 
two either false alleles or contaminations. The error rate in 
this study was low (see Bonin et al. 2004) and therefore un- 
likely to have notably influenced results from any analysis of 
the dataset. 

Descriptive genetics 

No statistically significant differences in genetic composition 
between the two cohorts sampled in 2006 and 2007 were ob- 
served in any of the sampling locations or between the tempo- 
rally spaced Faughan and Deele samples. Thus, all temporal 
samples representing all 13 sampling locations were pooled 
prior to subsequent analyses. 

The mean number of alleles per locus varied between 12.9 
in sample Owenreagh East and 18.4 in sample Faughan. The 
average expected heterozygosity (Table 3) was very simi- 
lar over all samples and ranged between 0.84 (Owenreagh 



South) and 0.88 (Fairy Water). Average allelic richness over 
loci (Table 3) varied from 1 1.99 (Owenreagh South) to 14.83 
(Faughan). While a few instances of significant (a = 0.05) 
genotypic linkage disequilibrium were detected in a few sam- 
ples, none of these were found to be relevant after Bonfer- 
roni correction for multiple tests. Similarly, only very few 
cases involving deviation from HWE were observed; and 
these involved not more than two loci per sample. None of 
these were found to remain significant following Bonferroni 
correction. 

Genetic differentiation, estimated as pairwise f st and Dest 
between sampling locations, ranged from 0.001 to 0.06 for 
Fsr and from 0.003 to 0.3 for Dest (Table 4). In general, Dest 
estimates were found to be in the order of three to four times 
higher than their Fst equivalents. All pairwise comparisons 
were found to be statistically significant with one exception 
involving samples from the Finn and Reelan. Global Fsr es- 
timate was 0.02, while the global Dest estimate (calculated as 
an approximation of the harmonic mean of Dest values over 
all loci using the following equation: Hmean ~ 1/ [(1/A) -I- 
var(D)(l/A)^] where A is the arithmetic mean and var(D) 
is the variance of Dest values over loci) was 0.12 (i.e., sixfold 
higher). Despite ongoing debate over the application Jost's D, 
and the merits of using _Fst for microsatellite data, we found 
a near-perfect correlation between both estimates (r = 0.96, 
P < 0.001) in the present dataset. Interesting, while we also 
found a strong positive correlation between Dc and Dest (>' = 
0.81, P < 0.0001), this was not as strong as between Fsr and 

Geographical population structure 

The NJ tree based on the Chord (Dc) genetic distance mea- 
sure (Fig. 2) confirmed the existence of several population 
groupings within the Foyle catchment, which were supported 
by moderate/high bootstrap support. These groups gener- 
ally adhere to the geographical origin of particular samples 
within the River Foyle system, thus samples from the west- 
ern, northeastern, and southeastern tributaries fall into three 
reasonably well-defined clusters. One of these geographical 
groups (western), identified by both NJ and STRUCTURE 
analyses, consists of geographically proximate samples from 
the western of the catchment. While the other two groups are 
characterized by geographically close samples (e.g., Faughan 
and Burn Dennet in the north and the Owenreagh South and 
Drumragh in the south) from the extreme ends of the catch- 
ment (i.e., south and north of the catchment), the analyses 
include other samples that are geographically inconsistent 
with the groups (e.g., eastern samples Camowen and Strule 
grouping with the north and the eastern samples of Glenelly 
and Owenreagh East grouping with the south). 

There was significant overall positive association between 
geographical distance and Fsr genetic distance among the 



364 



® 201 1 The Authors. Published by Blackwell Publishing Ltd. 



D. Ensing etal. 



Disentangling Evolutionary Signal from Population Stochasticity in Atlantic Salmon 



Table 3. Expected (He) and observed heterozygosity (Ho), number of alleles per locus (A/a), allelic richness (Ar), and inbreeding coefficient (fis). 
Samples showing significant deviation from Hardy-Weinberg equilibrium (HWE) are denoted with an asterisk (*). 
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Table 3. Continued 
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samples shown by the Mantel test (r = 0.293, P = 0.039; 

Table 5A). Not surprising, a similar positive association was 
also observed with Dest (r = 0.332, P = 0.009). This observed 
overall positive IBD pattern, however, was strongly influ- 
enced by the major regions (i.e., western, northeastern, and 
southeastern). Within each individual region, there was no 
evidence of IBD, vdth the exception of the western regional 
group. In the latter case, however, given small number of 
samples within the group, results should be interpreted with 
caution. To test the independent effects of both geographical 
distances and the major geographical regions, as identified 
by the previous analyses (e.g., STRUCTURE), in the patterns 
of population differentiation observed, partial Mantel tests 
were carried out (Table 5B). Genetic distances (Fsj and Dest) 
were found to be significantly correlated to both geographical 
distances and regions, with the latter variable having a more 
significant effect. Together, these analyses provide unambigu- 
ous evidence for IBD in the River Foyle among regions, but 
not within all regions per se. 

The results of the STRUCTURE analysis are in agreement 
with the Chord genetic distance clustering in supporting three 
major genetic groups within the Foyle catchment, which are 
generally defined by geography (i.e., western, northeastern, 
and southeastern tributaries). Interestingly, at the most basic 
level of structure, AK was found to be bimodal at IC = 3 
and K = 4. Close examination of the graphical output of the 
STRUCTURE analyses (Fig. 3) suggests some evidence for 
further partitioning within the northeastern group. Thus, 
the eastern samples from the Owenreagh East, Camowen, 
and Strule appear to be somewhat genetically distinct from 
samples from the Faughan, Deele, and the Dennet. Addi- 
tional STRUCTURE hierarchical analyses, however, failed to 
indentify additional structuring within this group. This was 
not the case for the southeastern cluster in which further 
substructuring was uncovered from the hierarchical STRUC- 
TURE approach. Thus, the Owenreagh, the Drumragh, and 
the GleneUy are clearly genetically distinct units (Fig. 3). The 
results from STRUCTURE and the Dq clustering were also 
confirmed by the Barrier analysis in separating the three ma- 
jor regional groupings within the Foyle System and also, to 
a large extent, the substructuring within the groups (results 
not shown). 



Geographical partition of genetic variance 

The results of the AMOVA analysis, carried taking in con- 
sideration (1) the main regions within the Foyle System (i.e., 
west, northeast, and southeast); (2) samples within regions; 
and (3) samples within the Foyle System, are summarized in 
Table 6. While 97.7% of the observed genetic variation is con- 
tained within samples, a significant (P < 0.001) percentage 
of the observed genetic variance is explained by differences 
among samples within (1.35%), and among regions (0.98%), 
thus confirming previous analyses. 

Discussion 

The overall levels of genetic variability in salmon of the Foyle 
River system (measured as heterozygosity, allelic richness, 
and mean number of alleles per locus) reported in this study 
are at the high end of the range reported for similar popula- 
tion genetic studies using microsatellite markers on Atlantic 
salmon (e.g.. Primmer et al. 2006; DiUane et al. 2007; Garant 
et al. 2000). This is possibly associated with the very large 
Atlantic salmon stock within the Foyle River system. 

Similar to previous studies examining within river system 
population structuring (e.g., DiUane et al. 2008), substan- 
tial levels of genetic differentiation were observed among the 
River Foyle tributary samples. Given lack of differentiation 
between the temporal samples examined in this study, it is 
reasonable to assume that the patterns of population genetic 
structuring observed within the Foyle River system are tem- 
porally stable. It is worth noting that one of the temporal 
samples was collected from a relatively small tributary in 
which the resident salmon population is believed to be small 
(i.e., average of 62 annual spawners over a 10-year period. 
Loughs Agency Annual Report 2005). 

The higher levels of genetic differentiation revealed by Dest 
relative to Fst are in hne with Jost (2008) , who argued that Fst 
(Gst) tends to underestimate true levels of population dif- 
ferentiation compared to Jost's D and related statistics (e.g., 
Hedrick 2005). Both methods, however, were in agreement 
with respect to the patterns of population structure observed 
within the Foyle system. The differences in the magnitude of 
the observed differences, however, could have implications 



366 



® 201 1 The Authors. Published by Blackwell Publishing Ltd. 



D. Ensing etal. 



Disentangling Evolutionary Signal from Population Stochasticity in Atlantic Salmon 



E 
U 



o o o o 
o o d o 



o o o o 



CO rsi 
o 
o 
d 



o o 
d d 



o o 
d d 



pn Lo 



d d 



tD l\ ^ O 



d 



o o o o 
d d d d 



o o o 

d d d 

— K r-v 

CM V- CM 

.~ rSi ^ 

S> S> S> 



CO 




fN 










Ln 






CO 


cri 






"X) 


Ol 


o 








O 






o 


o 






O 


o 


O 


o 


o 


o 


o 


o 




fN 


CO 


LD 


rsl 


Ln 


o 


m 


m 




^ 


rsl 




fN 












O 


o 


O 


O 


O 


o 


o 






o 


d 


d 


d 


d 


d 


d 


d 


1 








CO 


CO 




IN, 














rv-| 




r\j 








o 


o 


o 


O 


o 


O 






O 






o 


O 


o 


O 








LD 


o 


Ln 


00 


CO 


r~-- 




CO 


m 








o 


o 


o 




cn 


o 


O 


o 


o 


o 


o 


o 




rsj 


d 


d 


d 


d 


d 


d 


d 


1 












o 


fV^ 




























o 






o 












§ 


















o 




(N 


00 


o 








m 


CM 




(N 










rsj 


o 


o 


o 


O 


o 


o 




oi 


rsl 


d 


d 


d 


d 


d 


d 


1 


rsl 


cn 


o 


r\ 


CO 
















i\ 


















o 


O 














O 


o 


O 


O 


o 










O 




















fN 




fN 


o 






00 


fN 




















d 


d 


d 


d 


d 


1 


m 




cn 


o 


















Ln 




o 


rn 
















o 


c3 












o 


O 




o 














LTl 




<Ti 




<Ti 


Ln 


CO 


rsj 


m 


CM 


O 


O 




O 






LD 


o 


o 


o 


o 




cn 


uS 


(N 


d 


d 


d 


d 


d 


1 


m 


00 




LD 
























o 


















o 














o 


O 


o 






















no 


Gl 






CO 




m 


o 




cn 




d 


ro 


rsj 










iri 


d 




oi 


rsj 


d 


d 


o 


1 


Ol 


CO 






1^ 


















































































cn 


Ln 


m 


m 


m 


s 


8 








Ln 




LD 


o 










rM 






Ln 


cn 


d 


O 










cn 




m 










































































ID 






LT) 






'sT 


rM 


r-- 


s 




r-- 


fN 




UD 


fN 




cn 






CO 


fN 


d 


fN 




LT) 












cn 




5 


ro 


cn 



rsi 

^ Lo 

m rsj uo m 

o o o o 



CO LjD ro -nI- cn ^ rn 
^(-jib'^^^'^rsi 
r^orsju-i-^r-'CTio 
no'^'^oi'^^ro^ 



nn Ol CO 

tN "JD 

o o o 

o o o 



o o o I 



ro ro ^ CO 
r^ rsl ^ ^ 



rM cn LD LD 



_ _ rsj LD ^ 
Oir^-corNi'—^fYi 

LD^cdr-^'-^oirv-irorod 
LOcncoroLDLnLnco*DPO 



Ol r-. m 00 

Ln rM LD LD LD 



1— LDroroLncorsir-- 

Ln LO cn 

cri 00 cri 

cn Ol CO 



_ OJ i/i 

^ -3! S 5 5 a != 

i£ U Q O O Q iZ 



from a management viewpoint if conservation of population 
structure and genetic differentiation are the focus. For in- 
stance, on the basis of low Fsr estimates, populations that 
would look only marginally differentiated and perhaps not 
worthy of discrete population status and possibly associated 
specific conservation measures, would have to be classified as 
highly differentiated on the basis of D^st and therefore might 
justify a higher level of protection. 

It is important to note that, in many ways, this discussion 
arises because of a misunderstanding of Fsr statistics as it 
relates to microsatellite data, as opposed to low or moder- 
ately polymorphic allozyme data for which this statistic was 
particularly suited for. As clearly demonstrated by Hedrick 
(2005) and Jost (2008), standard Fsr estimates from highly 
polymorphic microsatellite data are often biased toward the 
lower end of the differentiation scale (i.e., 0 = nondifferen- 
tiated; to 1 = completely differentiated populations). Thus, 
there is an inherent methodological penalty associated with 
the use of this class of markers with this statistic. Many 
authors, however, stiU fail to acknowledge or recognize this 
important issue, with the unfortunate consequence that sig- 
nificant genetic differences are often dismissed as being in- 
consequential on the basis of "small" Fst estimates. Both 
Hedrick (2005) and Jost (2008) have suggested potential ad- 
justments to this statistic to account for the bias. This, how- 
ever, has not as yet been generally implemented in studies 
that have been reported in the literature. A potential rea- 
son for this is the fact that little is stiU known about the 
behavior of these "adjusted" Fsr estimates in relation to 
the evolutionary factors they are measuring, for example, 
gene flow, mutation, random genetic drift, and selection (but 
see Whitlock 2011). Thus, while traditional Fsr statistics can 
be readily interpreted in terms of the above, this is not nec- 
essarily the case for the "adjusted" Fsr statistics. Recendy, 
Whitlock (2011) suggested that the traditional Fst statistics 
perform better than the "adjusted" equivalent measures. In 
here, however, based upon evaluation of empirical data, we 
found no evidence supporting this argument beyond the re- 
ported differences in magnitude. Thus, further studies are 
StiU urgendy required to address this important issue. 

In instances where these "adjusted" measures of genetic 
differentiation have been implemented, results have been in- 
conclusive. For instance, Hoffman et al. (2009) did not find 
substantial evidence to support the better performance of 
Jost's Dest as compared to standard Fsr in detecting differ- 
entiation in Steller's sea lion (Eumetopias jubatus). We ar- 
gue, however, that levels of microsatellite polymorphisms 
observed in the Steller's sea lion (4-11 alleles over 13 loci) 
are comparatively low in relation to what is observed in other 
species (e.g., 13-37 alleles over seven loci in the present study 
for Atlantic salmon). Thus, the bias associated with high lev- 
els of genetic diversity as discussed by Hedrick (2005) and 
Jost (2008) was potentially less evident in Hoffman et al's. 
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id 



Owenreagh S. Qreen = Southeast. 



Table 5. (A) Mantel test results for correlation between genetic (FST and Dest) and geograpliic distances for botin the Foyle system and for tine major 
geographical groups identified during analysis; (B) Partial Mantel test results to investigate the independent regional and geographic effects on the 
genetic distance (FST and Dest). P-values derived from 10,000 iterations. Significant results (P < 0.05) given in bold. 

Group FsT Dest 



r-value P-value r-value P-value 

(A) 

All 13 2006-2007 sampling sites 0.293 0.039 0.332 0.009 

Northeast sampling sites only -0.142 0.692 -0.193 0.757 

Southeast sampling sites only -0.790 0.828 0.237 0.658 

West sampling sites only 0.729 <0.001 0.850 <0.001 

(B) 

Partial Mantel test 

Genetic and geographical distance controlling for regional effects 0.336 0.033 0.340 0.023 

Genetic distance and region controlling for geographical distance effects 0.431 <0.001 0.470 <0.001 



(2009) study in comparison to what was observed for At- 
lantic salmon in River Foyle. 

As conservation and management strategies of many 
salmonid fish species have been predominantly based on pop- 
ulation differentiation estimates of Pst> as recommended by 
Heller and Siegismund (2009), we argue that the use of ei- 
ther Hedrick's (2005) standardized genetic statistics or Jost's 
(2008)Dest should provide a more appropriate measurement 



of population differentiation. Indeed, we suggest that reviews 
should be carried out to reassess the levels and patterns of 
population divergence from studies solely relying on uncor- 
rected PsT- This should be a priority for valuable Atlantic 
salmon populations where the ability to accurately identify 
individual population units is critical. Furthermore, from a 
conservation perspective, it would be relevant to inform man- 
agers about the caveats of the interpretation of standard Fsr 
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Table 6. AMOVA among the Atlantic salmon samples from the River 
Foyle system. 
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Figure 3. STRUCTURE bar plots (aligned using CLUMPP) depicting pop- 
ulation structure based on two hierarchical rounds of analysis; K = 3 
(first round, left) and K = 3 (second round, right) where membership 
of each individual to the different clusters can be deducted from the 
proportion of each color in the bars. 



previously estimated and presented insofar as microsatellite 
data are concerned. 

Notwithstanding the issue discussed above, for a compar- 
ative perspective between this and other genetic studies on 
Irish Atlantic salmon populations using similar microsatellite 
markers, it is worth noting that the global uncorrected f sx 
values (f ST = 0.02) did not differ greatly. For instance, Dillane 
et al. (2008) reported global f st values of 0.024 for the Moy 
River in the west of Ireland. D^^t estimates, in the Foyle, both 
over all loci and population samples were on average three 



0) 
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P-value 


Among regions 
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<0.001 


Among samples within regions 
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<0.001 
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2157 
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to four times higher compared to global f st values. Again, 
using the argument above, if this study had been limited to 
only using Fsj as a measure of genetic differentiation, the 
vast majority of differentiation between samples could have 
been perceived as very small. Given the argument above, the 
important point is whether the statistic is significant or not, 
while the value itself is largely irrelevant. 

Given observed differences in allelic frequency distribution 
(exact tests), the majority of the samples examined from the 
Foyle represent more or less independent temporally stable 
population units (i.e., with limited levels of gene flow). The 
only exception (i.e., no genetic differences) involved the geo- 
graphically close sampling locations Finn and Reelan (west- 
ern region, Fig. 1). Given that these samples were taken from 
within the same tributary river, the lack of genetic differenti- 
ation between them is not entirely surprising. 

Despite genetic differences observed among individual 
samples, our data suggest the existence of three major pop- 
ulation groupings (West, Southeast, and Northeast) in the 
Foyle catchment, the overall pattern of which generally ac- 
cords with the geographical and hydrological organization 
of the catchment, and which is consistent with a member- 
vagrant evolutionary model. However, while the membership 
of samples belonging to the West group is geographically well 
defined, some inconsistencies were observed in the member- 
ship of the Northeastern and Southeastern groups. We argue 
below that these discrepancies can be explained by anthro- 
pogenically mediated disturbances. 

The western tributaries comprise the largest rivers in the 
Foyle catchment with correspondingly large population sizes. 
In comparison to the other areas of the catchment, the west- 
ern rivers are considered to be reasonably pristine from an 
environmental perspective. This area of the catchment is also 
well known for the presence of multi-sea winter fish, a life- 
history trait that is uncommon in the rest of the catchment. 
Thus, the genetic uniqueness of samples from this area is 
likely to reflect an independent evolutionary history. Overall, 
there was no strong evidence for IBD within groups, which 
suggests high levels of conductivity and therefore is in agree- 
ment with the concept of a metapopulation evolutionary 
model. 

Within the Northeast group, the most northern tributaries 
(i.e., Faughan, Deele, and Burn Dennet), similar to those in 
the west of the catchment, are also large, with good salmon 
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habitat and supporting large populations. Indeed, available 
demographic data for the Foyle system (authors' unpublished 
data) suggest that the tributaries from the north and from 
the west of the catchment are responsible for the most of 
the current salmon productivity in the Foyle (i.e., ~80%). 
In contrast, however, the population sizes of the eastern 
(GleneUy, Owenreagh East, Camowen, and the Strule) and 
southeastern (Owenreagh South and Drumragh) tributaries, 
despite their large habitat area (between 25% and 35% of the 
catchment) are relatively narrow and small (authors' unpub- 
lished data). This part of the Foyle catchment is character- 
ized by intensive farming and it is heavily impacted by arterial 
drainage schemes resulting in the disturbance of natural habi- 
tats, which severely Hmit productivity. Therefore, the genetic 
pattern characteristic of samples belonging to this region re- 
flects random stochastic processes (e.g., genetic drift acting 
on small gene pools), associated with anthropogenetic im- 
pacts, such a loss of spawning habitats or poor water quality, 
as the main factors shaping contemporary population genetic 
structure in these regional groupings. 

This hypothesis is supported by the lack of evidence for 
substructuring as provided by both NJ and the STRUCTURE 
analyses. Alternatively, the observed similar genetic makeup 
of samples from both regions could reflect past historical 
gene flow prior to the development of the extensive arte- 
rial drainage scheme in the eastern and southern regions of 
the catchment. At the moment, the data are insufficient to 
disentangle both hypotheses. 

While we observed strong and statistically significant pos- 
itive correlation between geographical distance and genetic 
distance (IBD) within the River Foyle Atlantic salmon stock, 
this trend was Hnked to the three indentified genetic regional 
groupings rather than individual samples. Interestingly, IBD 
within a single river system is relatively uncommon in At- 
lantic salmon. Evidence for IBD is more frequently reported 
from studies on Atlantic salmon population genetic structure 
at larger spatial scales often involving multiple river catch- 
ments (e.g., McConnell et al. 1997; King et al. 2001; and 
Dillaneetal. 2007). 

The Atlantic salmon stock of the River Foyle most likely 
exists as several distinct, long-term temporally stable, pop- 
ulations within well-defined geographical regions, in which 
migration occurs more frequently than between populations 
from different regions. It is clear from our findings that the 
genetic structure of the River Foyle Atlantic salmon stock 
is complex. Thus, it cannot be simply defined within a sin- 
gle evolutionary model. If considering the whole system, the 
population structure of salmon within the Foyle catchment 
is better explained under the member-vagrant evolutionary 
model as confirmed by the clear presence of IBD at this level. 
However, within each of the major genetic regions identified, 
the metapopulation model best explains the observed struc- 
turing, as there is no evidence for IBD at this level. We argue 



that these models are not mutually exclusive and can op- 
erate simultaneously on different geographical scales within 
a single river system. Thus, any attempts to identif)^ a best 
model fit would be artificial and not reflective of the complex 
evolutionary history commonly associated with the species. 

Interestingly, Garant et al. (2000) has previously argued 
that the member- vagrant evolutionary model was too rigid to 
explain the population genetic structure of Sainte-Marguerite 
River salmon populations. In their study, the authors found 
difficult to accommodate extinction-colonization events, 
lack of temporal genetic stability, and absence of an IBD with a 
pure member- vagrant hypothesis, as these processes are more 
associated with the metapopulation model. Similarly, while 
Primmer et al. (2006) concluded that the member-vagrant 
evolutionary model best explained the genetic structure ob- 
served in the Varzuga River, they have also reported possible 
extinction-colonization events prompting a consideration of 
metapopulation processes acting at some sites. 

From a conservation and management viewpoint, focus 
should be orientated toward identifying and preserving ap- 
propriate evolutionary entities. In the Foyle, using genetic 
evidence and an environment quality indicator (EQI), it 
is evident that conservation efforts should be targeted to- 
ward southern and eastern tributary catchments. This could 
be done by restoring the integrity of habitats prior to the 
development of intensive farming activities and associated 
drainage schemes. The data also suggest that the manage- 
ment of salmon in the Foyle should accommodate the three 
regional groupings identified in this study. Future stud- 
ies should continue to explore the interplay between the 
evolutionary models determining genetic structure and 
to what extent the findings reported in here apply 
elsewhere. 

Data Archiving 

Genotypes and all other raw data used in this study are stored 
on an electronic database at AFBINI. At request available to 
others via the corresponding author. 
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